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,_H . We consider a special class of the tetrad theory of gravitation which can be considered 

£Nj ■ as a viable alternative gravitational theories. We investigate cosmological models based on 

those theories by examining the possibility of fitting the recent astronomical measurement 
of supernova la magnitude versus shift. Our investigations result in a reasonable fit for the 
supernova data without introducing a cosmological constant. Thus, cosmological models 
based on tetradic theory of gravitation can provide alternatives to dark energy models. 
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1 Introduction 

The present observation of the distant supernovae type la indicates that the universe is presently 
acceleratingfl]— [4]. The cosmic acceleration is attributed to the presence of unknown form of 
energy violating the strong energy conditions px + 3px > where px and px are energy density 
and pressure of dark energy respectively. Different candidates for dark energy are attempted to 
yield accelerating cosmologies at late time[5]-[14]. The cosmological constant A and phantom 
fields [15]— [17] violating a weak energy conditions px + Px > are most popular ones. 

Recently, different attempts [18]-[25] have been carried out to modify gravity to yield accel- 
erating cosmologies at late times. In this paper, we exploit the possibility of modifying gravity 
based on absolute parallelism spaces. 

The notion of absolute parallelism was first introduced in physics by Einstein [26] trying to 
unify gravitation and electromagnetism into 16 degrees of freedom of the tetrads. His attempt 
failed, however, because there was no Schwarzchild solution in his field equations. 

The interest in the tetrad theory as a purely gravitational theory was revived by M0ller[27] 
who showed that a more satisfactory treatment of energy momentum complex than that of 
general relativity can be achieved. In his first attempt in finding Lagrangians Mpller's was 
restricted by the assumption that the equations determining the metric tensor should coincide 
with the Einstein equations. After then, he [28] abandoned this assumption and relooked for a 
wider class of Lagrangians, allowing for possible deviation from the Einstein equations in the 
case of strong gravitational fields. Mpller's theory was generalized into scalar tetradic theory by 
Saez[29]. Meyer showed that Moller's theory is a special case of the Poincare gauge theory [30]. 



Quite independently, Hayashi and Nakano[31] formulated the tetrad theory of gravitation 
as a gauge theory of the space-time translation group. Hayashi and Shirafuji[32] studied the 
geometrical and observational basis of the tetrad theory, assuming that the Lagrangian be given 
by a quadratic form of torsion tensor. If the invariance under the parity operations is assumed, 
the most general Lagrangian consists of three terms with three unknown parameters to be fixed 
by experiment. Two of the three parameters were determined by comparison with solar-system 
experiments, while only an upper bound has been estimated for the third [32]-[33]. 

The numerical values of the two parameters found were very small, consistent with a value of 
zero. If these two parameters are exactly equal to zero, the theory reduces to the one proposed 
by Hayashi and Nakano and M0ller which we shall here refer to as the HNM theory. This theory 
differs from general relativity only when the torsion tensor has a nonvanishing axial vector part. 

Many applications of HNM theory have been done during the years. These include cos- 
mological applications [34], investigating gravitational radiation [35]— [36] and energy momentum 
complex [37]-[38], and finding a general solution with spherical symmetry [39], and solution with 
axial symmetry in [40]. 

The rest of the paper is organized as follows: In section 2 we review HNM tetrad theory of 
gravitation. In section 3 we present the tetrad field satisfying the requirement of homogeneity 
and isotropy. In section 4 we discuss briefly the basic of the tetrad cosmology and derive 
the relation between the luminosity distance and redshift. In section 5 we present numerical 
investigations of the tetrad cosmological model to fit the recent observational supernovae data. 
The obtained results are compared with standard cosmology (with or without cosmological 
constant). Finally, in the last section we give our discussion. 



2 HNM Tetrad Theory of Gravitation 

In this paper we follow M0ller construction [28] of the tetrad theory of gravitation based on the 
Weitzenb0ck space-time. In this theory the field variables are the 16 tetrad components e^, 
from which the metric is constructed as 

<T := V ij ei%\ (1) 

where rf 3 is the Lorentz metric tensor taken as diag(— 1, 1, 1, 1). The Latin indices . . .) refer 
to vector numbers and Greek indices (/x, v . . .) to vector components, and all of them run from 
to 3. We restrict indices a, b, . . . and a, (3, . . . (beginning of Latin and Greek alphabetic) for 
spatial components. 

An invariant Lagrangian C is constructed from g tlv and 7 ^ p , where 7 ^ p is the contorsion 
tensor given by: 

Ipvp '■= T\ 3 ej^ Gju-p, (2) 

where the semicolon denotes covariant differentiation using the Christoffel symbols. The most 
general Lagrangian density invariant under the parity operation can be constructed as a linear 
combination of the following expressions: 

L« := Q^, : = 7/w7 ^, : = 7 _ 7 ^, (3) 

where 3> M is the basic vector defined by 

== 7V ( 4 ) 

These expressions in eq. (3) are homogeneous quadratic functions in the first order deriva- 
tives of the tetrad field components. 



M0ller considered the simplest case, in which the Lagrangian C is a linear combination of 
the quantities i.e., the Lagrangian density is given by 

^M0ller := ("S) 1/2 + a 2 L< 2 > + « 3 ^ (3) ), (5) 

where 

g := det(g^). (6) 

For this choice, the constants on had been chosen such that this theory gives the same results 
as GR in the linear approximation of weak fields. According to his calculations, one can easily 
see that if one chooses 

Qi = —1, «2 = A, 03 = 1 — 2A, (7) 

with A equals to a free dimensionless parameter of order unity, the theory will be in agreement 
with GR to the first order of approximation. The same identification of the parameters was 
obtained by Hayashi and Nakano[31] 

M0ller applied the action principle to the Lagrangian density eq. (5) and obtained the field 
equations in the form [28] 

= 0, (8) 
where the Einstein tensor G^ v is defined by 

G^u := R^u — - g^ u R- (9) 
Here and are given by 



Hfj,u '■— A 



Fpv '■— A 



'laP^u + Jaf3^» a(B + Jaf3»^ a(B + 9^ (lafial^ ~ \la^ C ) } (10) 

and 

%,* ~ ~ $>* (tV " 7%) + 7^%] • (11) 



The term H^ u by which equations eq. (8) deviate from Einstein's field equations increases with 
A, which can be taken of order unity without destroying the first order agreement with Einstein's 
theory in case of weak fields. 

M0ller assumed that the energy-momentum tensor of matter fields is symmetric . In the 
Hayashi-Nakano theory, however, the energy-momentum tensor of spin-1/2 fundamental parti- 
cles has a nonvanishing antisymetric part arising from the effects due to the intrinsic spin. 



3 Tetrad fields for applications to cosmology 

The tetrad fields satisfying the symmetry requirements of homogeneity and isotropy have been 
given by Robbertson[41]. It was found that there are two possible teterads, which in Cartesian 
coordinate can be written in the form 

e ° = 1, e a ° = 0, e a = 0, 

R(t)e a a = 5 a a h^ + lx a x a ±khe a ^x^ (12) 



e ° = ^, e a ° = ±^-^x a , R(t)e a = ±(-k)^x a , 

R(t)e a a = 5 a a h + -^x a x a , (13) 

where the constant k takes the values +1,-1, or zero, while h± and r 2 are defined by 

h± = l±jr 2 , r 2 = x 2 + y 2 + z 2 . (14) 

Here, e aa fi = ±1 when (cm/3) is an even or odd permutation of (123) and otherwise. 
Both of these two tetrads through eq. (1) lead to Rob erston- Walker metric given by 

ds 2 = -dt 2 + R{tf h^ 2 {dx 2 + dy 2 + dz 2 ), (15) 

Explicit calculations based on eqs. (2)— (4) and the Christoffel symbols of the metric defined by 
eq. (1) result in the following nonvanishing components of the tensor 7 and $ for the first tetrad 
in eq. (12). 

7 ° afe = - 5 ab R R h+ 2 , 7 « 6 = - 5 ab |, 

lie = - tabc ^ $0 = - 3|, $ a = 0, 

(16) 

where R =§ §- 

Concerning the second tetrad, because the physically relevant second rank tensors derived 
from it, could be complex as pointed out in [42]. Thus, it will be dropped out in the present 
work. 



4 The basic equations of cosmology 

The skew symmetric part of the field equations in eq. (8) is satisfied identically by the tetrad de- 
fined by eq. (12), which can be verified by explicit calculations. Assuming the energy-momentum 
tensor for a perfect fluid, then the symmetric part of the field equations in eq. (8) reduces to 

(2R/R) + (R/R) 2 + D/R 2 = -p, (17) 

3(R/Rf + ZD/R 2 =p, (18) 

where D = k (1 — 3 A), p is the pressure, and p is the density of energy associated with the fluids. 
Through out this work we employ units in which G = 1/8tv and c = 1, and in this section we 
use the notations and equations (with suitable modification) presented in [43] 

As in the standard cosmology, the Hubble parameter H is defined by H = R/R. Then 
critical density is given as 

Pc = 3 (i?/i?) 2 = 3H 2 (19) 

Many useful parameters can be defined, among them the density fi and the deceleration param- 
eters defined respectively as: 

n = p/p c , q = -RR/R 2 , (20) 

In this manner, eq. (17) and eq. (18) take the form 



q = n/2 + P /2H 2 , n + n D = i, 



(21) 



where Qd = —D/R 2 H 2 . It is worthy to mention that the metric in eq. (15) and the field 
equations eqs. (17)-(18) involve two different constants, namely k and D respectively. Since A 
is completely free, then the two constants are independent. Thus, there is no unique relation 
between the geometry of the universe and its fate, in contrast to the case of general relativity 
(without cosmological constant), where such a relation exists. 

For matter dominated universe eq. (18) can be written in the form 

PoRq 



R z 



R 



-D, 



(22) 



here D can take any real value in contrast to the case of general relativity where k is restricted 
to ±1,0 values. 

Let us define the parameters 

1 R 



H t 



1 + z Rq ' 

where z is the red shift. In terms of these parameters eq. (22) can be transformed into 



da^ 2 
dr , 



l + ^o 



- 1 



which can be converted, in terms of z into 



dz 
~dr~ 



(i + zf (i + n m0 z)2. 



(23) 



(24) 



(25) 



The time lapse between the present time and the emission time of light which suffers redshift z\ 
can be obtained from eq. (25) leading to 

rzi 



to — t\ 



rzl I 

Jo + , 



(i + n m0 z) 2 dz. 



Light signal moves along null geodesic whose equation is 
(1 - kr 2 Y 



dr 
dt 



R 



and using eqs. (23)-(25) we obtain 



dr 



[l - kr 



2\V2 



[i + z) = {i + n m0 z) 



Rq 



RoH {1 + z) 



dz. 



(26) 



(27) 



(28) 



Integrating eq. (28) for a light traveling through the universe from r\ at time t\ , and reaching 
at r\ = at the present time to we obtain 



rri 

Jo (T 



dr 



where 



rT 

Jo 



kr 



dr 



2\V2 



zi 



(1 + 17 m0 z) 

R H (l + z) 



dz 



(29) 



(1 - kr 2 ) 



2\V2 



sinh^ 1 (|fc|) 



(1*1) 
l 



T72 



1/2 



(1*1) 

n 



Tj2 sin 1 /-j 



.1/2 



r\ k 
k 



k 



1 





(30) 



From Qdo = —D/RqHq, we can express \k\ as 
\Vdo\ R 2 H 2 



\k\ 



(31) 



where 6a = 1 — 3 A. Using eq. (29), for the case 6a > and k = — 1, we obtain 



1/2 



> I ^DO I / RqHq 

then solving for r\ one gets 



sinh 1 



\9. 



DO 



1/2 



Ro H r\ 



ji + n m0 z) i 

o RoH (l + z) 



dz 



n 



1/2 



1 



■ sinh 



\n 



DO 



1/2 



2 i (i + n m0 z) 



\^Do\J RqHq " " \ V |^a| / ./0 

The luminosity distance can be shown to be 
d L = Rl n/Ry, 
and with the help of eq. (33) we get 



(1 



- dz 



/ M \" 2 (i + z) . 
d ' = (,^J — smh 



10 



DO 



1/2 



(i + rt m0 u) 

(l + n) 



■ du 



(32) 



(33) 



(34) 



(35) 



5 Testing the model against the Supernovae data 

The supernovae of type la (SNe la) serve as excellent cosmological standard candles. The 
apparent magnitude of a "standard candle" is related to its luminosity distance di through 



m(z) = M + 5 log 



10 



d L 
Mpc 



+ 25, 



(36) 



where M is the absolute magnitude and is assumed to be constant for a standard candle like 
Sne la. The apparent magnitude also can be expressed in terms of the dimensionless luminosity 
distance T>l(z) as 



m (z) = M + 5 \og 1Q V L (z), 



with 



and 



V L {z) = ^d L 



M = M + 5 log 10 ) + 25 



(37) 
(38) 



\lMpc 
M - 5 log 10 h + 42.38. 



(39) 



For our present analysis we use "gold" sample compiled by Reiss et al.[4]. The sample 
consists of 157 data points which in terms of distance modulus are 



^obs = m ( z ) ~ M > 

= 5 log 10 V L {z) - 5 log 10 h + 42.38. 

The best fit model to the observation is obtained by using x 2 statistics, i.e. 



X 



157 

E 

i=i 



^th ^obs 



-i 2 



(40) 



(41) 



where fj,^ is the predicted distance modulus for a supernova at redshfit z and Oi is the dispersion 
of the mesaured distance modulus due to intrinsic and observational uncertainties in SNe la peak 
luminosity. In our model we obtain the best fit for O m and 6a, from the minimization of x 2 by 
scanning the whole relevant parameter space, the values we obtain are O m = 0.3 and 6a = 0.26 
for k = — 1 and h = 65 km s _1 Mpc _1 . The other values of k don't lead to a best fit. In 
fig. 1 we illustrate the variation of x 2 versus 6a for fixed Q m = 0.3 and k = — 1. For a one 
parameter fit (6a), the Sne la data provides the following ranges: 0.234 < 6a < 0.288 at 69% 
confidence level (CL) and 0.207 < b\ < 0.314 at 95% CL. For a quantitative comparison we use 
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Figure 1: Variation of x 2 with 6a- 





Model 




x 2 




= 0.30, ft A 


= 0.70 


178 




= 0.30, Q A 


= 0.00 


273 




= 0.30, 6 A 


= 0.26 


212 



Table 1: Values of x 2 

the x 2 fit- in table 5 we present the x 2 values corresponding to best fit for the general relativity 
with and without cosmological constant [4], and present work. The x 2 value for the present 
work seems to give a reasonable fit of the Sne la data without the necessity of introducing the 
cosmological constant. In Fig. 2, we also make a comparison between these theories using the 
effective luminosity and the redshift. Qualitatively, as can be seen from Fig. 2 the tetrad theory 
fits the data better than general relativity without cosmological constant and it is almost close 
to the general relativity with cosmological constant for the current available data. 

6 Discussion 

In this paper we consider a special class of the tetrad theory of gravitation as a viable alternative 
gravitational theories. We have made use of recent measurements of supernovae la to compare 
the tetrad theory with general relativity with and without cosmological constant. As has been 
shown in the previous section the tetrad theory leads to a reasonable fit of the supernovae data 
without introducing cosmological constant. However, to make a best of the data we find that the 
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Figure 2: Effective magnitude versus redshift z. 



corresponding universe is open and non flat. This may implement a new scenario for inflationary 
open universe. Such a scenario based on general relativity has been previously discussed by [44]. 

The density parameter Qm = 0.30 of our work is the same as that of general relativity with 
f&A = 0.70, however in our case the cosmological constant is already absent. Concerning the age 
of the universe, our model gives the same value as that of general relativity (Qm = 0.30, Q\ = 0) 
which is OSIHq 1 . In the presence of cosmological constant (£Im = 0.30, Qa = 0.70), the age of 
the universe turns out to be little more around 0.96HQ 1 . 
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